BLACK HOLE SOLUTIONS WITH DILATONIC HAIR IN HIGHER CURVATURE 
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A new numerical integration method for examining a black hole structure was realized. Black 
hole solutions with dilatonic hair of 4D low energy effective SuperString Theory action with Gauss- 
Bonnet quadratic curvature contribution were studied, using this method, inside and outside the 
event horizon. Thermodynamical properties of this solution were also studied. 
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I. INTRODUCTION 



ly-^ , During last years a great interest for an investigation of low energy effective SuperString Theory action in four 
dimensions arised. Some researches [|l|-|4| found that the well-known solutions (such as Schwarzshild one and ets.) 
were modified by the higher order curvature corrections. They also showed that the "no-hair" theorem can not been 
applied to those modified configurations. 

The problem is to find black-hole solutions of 4D low energy String effective action with the second order curvature 
corrections. For simplicity, as a rule, researchers j^J^] consider the bosonic part of the gravitational action consisting 
of dilaton, graviton and Gauss-Bonnet (GB) terms taken in the following form: 



(1) 



where R is a scalar curvature; <f> is a dilaton field; m p i is the Plank mass; A is the string coupling parameter. The 
later describes GB contribution (Sgb = RijkiR^ kl — 4i?y i? y + R 2 ) to the action ([!]). Such configurations were partly 
studied g] by the perturbative analysis O(X) outside the event horizon when event horizon radius 3> m p i. Authors 
showed the black hole solution to be real and it provides the non-trivial dilatonic hair. P.Kanti et all || obtained the 
i similar solution using the non-perturbative numerical method outside the horizon. However, a question on a solution 
behavior inside the event horizon is still open, ft is generally believed || that in the regions where the spacetime 
curvature is just small a classical solution gives the main contribution to the global structure of the spacetime. 
Quantum corrections may drastically modify the spacetime properties in the case of large enough curvature. The 
purpose of this paper is to study the complete black-hole solution. 

Structure of our paper is the following: in Section 2 an analytical investigation of the action (^) is described; Section 
3 deals with the numerical results; Section 4 is devoted to description of thermodynamical properties of the solution; 
Section 5 contains a discussion and conclusions. In Appendix one can find a brief description of our integration 
method. 



II. ANALYTICAL INVESTIGATION OF THE GAUSS-BONNET ACTION 



The aim is to find static, asymptotically flat, spherically symmetric black-hole-like solutions. In this case the most 
convenient choice of metric is the following: 

2 

ds 2 = Adt 2 ~^dr 2 - f 2 (d9 2 + sin 2 Bdip 2 ), (2) 

where functions A, a and / depend only on radial coordinate r. Therefore, the scalar curvature R and the GB term 
Sgb have the following forms: 
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Integrating i) by parts the GB term and R term and ii) over the angle variables in ([!]) , one can rewrite the action in 
somewhat more convenient form (for the present analysis the boundary term is not relevant and is ignored): 



5' 



dtdr 



m%-{A'f'f + A(/') 2 + a 2 - A/ V) 2 ) + 4 e - 2 *A0'( — * ' -) 



Llf AA'{ff A 
a 



<7° 



Corresponding field equations in the curvature gauge (/(r) = r) are the following: 

m 2 pl a 2 [-ra' + ar 2 (<j>') 2 ] + 4e~ a *A(r(A - <r 2 )[0" - 2V) 2 ] + Ae~ 2 ^^Xa'(a 2 - 3A) = 0, 



(5) 
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- A'r - A] + 4e- 2 V AA'(cr 2 - 3A) = 0, 



m 2 pl a 2 [A"ra - AW + 2A'a - 2Aa' + 2Ara{(j)') 2 ] 

+4e- 2 *A2t7AA'[0" - 2(f) 2 ] + 4e- 2 VA2[(A')V + AAV - 3 AAV] = 0, 

-2m 2 i o- 2 [AV 2 o-0' + 2Ar(f)'a - Ar 2 <j)'a' + Ar 2 <j)"a] 
+4e- 2< ^A[(A') 2 CT + AAV - 3 AAV - AV 3 + AVer 2 ] = 0. 
The equation (Q) can be solved analytically and it permits to get the following expression for a: 
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In a formal mathematical way it is necessary to put the formula (|10|) to equations (|6j), (J8|), (|9|) and to work with two 
unknown functions A and <f>, but in the given situation it is more convenient to work with the equations (^]), (|^) 
and three unknown functions A, a and <fr, verifying (^) at every numerical integration step and writing asymptotic 
expressions for all four equations. We choose the later strategy. 

Let us suppose the solutions behaviors near the event horizon to have the following forms: 



A = dix + d 2 x 2 + 0(x 2 ), 
a = sq + six + 0(x), 
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(12) 
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where x = r — rfc,<S 1. Putting formula (jO)-([T^) in the equations (^)-(j^), we obtain the following equation and 
relations between the expansion coefficients [so, Si and 0o are free parameters): 



where: 



Zid x + z 2 dx + z 3 — 0, 



0/l\2i2 o \ 2 j.2 2 442 

zi = 24A <p Q r h , z 2 = 8A </> s - m pl r h s , 



(14) 



4 3 4 

z 3 = m pl r h s ; 



and other parameters are: 



2r h sl 



){rud\ + disl - rhsosidi], 



rn 



pi 



■[r h di - si], 



4\diq>o 

Q , 1 7 [8A(/) 2 ^o + 4Asi0i0 o + rls^lmlf - r h sim pl }. 



(15) 

(16) 
(17) 



2 



The solution of quadratic equation (|14|) exists if a value of the free parameter <f>o satisfies the following condition: 

x2 > 7(1 + VW5)m* pl 4 



2A 2 



(18) 



This is a natural restriction to the solution existence of the equation system @-(|^). 

We are interesting in asymptotically flat solutions. Therefore, its infinity behavior is to be as follows: 
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where ADM mass M and dilaton charge D can be calculated ^from the equations (0), (|S|) by the following formulae: 
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The dependencies of ADM mass M and dilaton charge D versus the event horizon radius are shown in Tab.l. 



III. NUMERICAL RESULTS 



For integrating inside the event horizon a method based on integrating over an additional parameter was used. It 
is briefly described in Appendix. Here we present the main results. 

The integration was done from the event horizon to the infinity, then, using data obtained, from the infinity to the 
horizon and, further, inside the horizon. The results of our calculation are shown in Fig.l. It presents the dependence 
of metric functions A and a and the dilaton function (f> versus r. Parameter Th was allowed to change within wide 
range from 3.5 to 100.0 Plank unit values. Fig. 2 represents 3D plot of the dependence of the metric function A(r) 
and dilaton function exp(— 2<p(r)) versus r%. As one can see from Fig.l and Fig. 2 the behaviors of A, a and cj) outside 
the horizon have the usual forms, as obtained in and look like the standard Schwarzshild solution (ADM mass 

values M (r/, ~ 2M) illustrate this, see Table 1). Under the horizon r < solution exist only till the value r = r s . 
Another solution branch begins from the value r s , but it exists only till the "singular" horizon r x . The asymptotic 
behavior of both solution branches near the position r s can be described by the following formulae, using the only 
smooth function a as an independent variable: 



A = d s + d 2 x 2 + 0(x 2 ); 
r = r s + r 2 x 2 + 0(x 2 ); 
exp(-20) = 0.(1 - 2f 2 x 2 ) + 0(x 2 ), 



(22) 



where x — a — a s <^ 1. Free parameters are the following: er s , 
be calculated from the following three equations: 

d2 = h, 
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The equation ( p4[ ) represents the asymptotic form of the equation (^) . The formula (|2^) represents the asymptotic 
form of the equation (||) and the rest of equations (§|) , (||) reduce into an identity. The equation (p5|) is the consequence 
of the system (|3^) (see Appendix) because according to the existence theorem (see Appendix) the system ( |35| ) has a 
single solution only in the case of its main discriminant to be not equal to zero. In the case the zero main discriminant 
in some point of the solution trajectory, the uniqueness of the solution ( |35| ) will be violated. So, equation ( p5| ) is just 
that condition in the r s , where two solutions exist. "Curvature invariant" RijkiK 1 ^ 1 in our metric parametrization is 
equal to: 



R ljk iR l3kl 
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According to the classification from r s represents pure scalar singularity. It is necessary to note, as we tested, 
that the weak energy condition, dominant one and strong one are realized at the value r — r s |6| . 

The most interesting result is the existence of the second solution branch within the "singular" horizon r x . The 
position of r x is situated inside the main horizon r^. The asymptotic behavior of A, a and </> near r x are the following: 

A = dix + d 2 x% + 0(x%), a = (j + aix? + O(x^), = <f> x + 0ix + fax^ + 0(x%), 



where x = r x — r <C 1, and free parameters are cto, <ji, 
from the following conditions: 
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"Curvature invariant" RijkiR^ kl ~ - + O(x^)) — > oo and this means that r x represents the "singular horizon". The 
distance between r x and is rather long for the big values of and decreases with decreasing rv In the limit point, 
defined by the equation (|l4|), all points pour together = r s — r x and the solutions inside do not exist. This can 
be seen from Fig. 2. 

IV. THERMODYNAMIC AL PROPERTIES OF GAUSS-BONNET SOLUTION 

If one works in Loretzian spacetime, he can say nothing about thermodynamical properties of a solution obtained. 
Using the Euclidean version of the metric: 



ds 2 = Adr 2 + —dr 2 + r 2 dfl, 
A 

where r is a periodic coordinate which range is from to tt , one can easily write the inverse temperature |p"i"|| : 

" d 
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^From equations (|ll|)-(|l3|) it is possible to find that (3 = 47r(so/di). 

In the Euclidean frame the full action Se after integrating over periodic time r takes the following form: 



Se = \p / dr 
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The boundary term obtained as it was done in M has the following form: 



2 —\ + Ae- 2 *\ 



(A A 1 



(29) 



4 



For a non-extreme black-hole the boundaries of the spacetime manifold are set by the extreme values of the radius 
coordinate. They are r = and r = oo. After removal of the conical singularity, the spacetime has the only boundary 
at r — oo, because the black hole horizon is not a spacetime border Q. So, we have to take into account only value 
of boundary term K. It takes the following form: 



K-K = (K- K ) rh = i/? 



'pi 



- 2r h ) - 40 O A— 
s s 



(30) 



in terms of the asymptotic solution on the event horizon (|T4|). Kq is the boundary contribution of the flat spacetime. 
In this case the Euclidean action Ie takes the following form: 



drS E -{K~ K ) r 



(31) 



where Se is defined by (E9). 

The physical entropy can be found by the following expression || : 



S(r h ) = [3^~I E = f3(r h ) 



dl(r h )\fdf3(r h ) 



dr h J \ dr h 



-I E (r h ). 



(32) 



Plots for the GB and Schwarzshild entropy are shown in Fig. 3. 



V. DISCUSSION AND CONCLUSIONS 



In this paper there are suggested the black hole solutions with non-trivial dilaton "hair" of low energy effective 
SuperString Theory with second order curvature corrections which are obtained independently from Kanti et all || 
and were found outside and inside the event horizon by the numerical method described in Appendix. The solutions 
are characterized by ADM mass M, dilaton charge D and asymptotic dilaton value </>oo. They are stable under 
the fluctuations of initial conditions. So, we can arrive to a conclusion that the SuperString Theory provides such 
situations by itself as can not be covered by the "no-hair" theorem |g. As these solutions have the non-perturbative 
nature, as they have no limits with the perturbative parameter values. 

The most interesting result of our work is the existence of r s -singularity inside the black hole after adding the GB 
term to the action. This singularity has the topology S 2 x R 1 , i.e. it is an infinite (in time direction) "tube" of radius 
r s . Similar "tube" in the Schwarzshild metric with an additional condition of R abcd R a bcd fmiteness was discussed by 
V.Frolov et al. There are two solutions on this "tube". The asymptotically flat solution, which is the main one, 
starts from r s and continues without limits up to the infinity. In the case of r% to be quite large and 3> r s it 
is possible to suppose r s to be approximately vanished, then the main solution looks like the standard Schwarzshild 
one with a constant dilatonic field (see Fig. 2) which agrees with results of Mignemi and Kanti ||. The additional 
solution branch provides the existence of the "singular" inner horizon with R^j^R? 3 — * oo. Some solutions exist 
inside "tube" r s , but they are unstable under the initial condition fluctuations, and we can not distinguish which 
branch — the main one or the additional one — they correspond to. One can suppose that this r s singularity is the 
consequence of the bad metric parametrization choice. In order to prove (or not) this artifact, the solution in another 
metric was studied. This is the metric of such class that was suggested by P.Breitenlohner et al. pl| and has the 
following form: 

ds 2 = f3 2 (r)dt 2 - dr 2 ~r 2 (T)d(l 2 . (33) 

Topology S 2 x R 1 inside the black hole horizon also exists in this metric. So the presence of r s "tube" is not the 
consequence of the metric parametrization choice. The question on the more detail structure of the "tube" r s is the 
subject of the next study. 

The Gauss-Bonnet term influence to the solution behavior was also considered. When the horizon radius r/j is quite 
large or A becomes quite small the main contribution to the action comes from the Einstein part (r/j 3> r s r s ~ 
in comparison with r^). While decreasing r/, (increasing A), the GB influence increases, the dilaton charge absolute 
value becomes larger, the distance between r s and becomes smaller. Once the limit position is achieved where 
fh = r s = fx = fhxs- This is a minimal point which provides the solution like the black hole type. In the point r) lxs 
the solution (main branch) exists only outside the event horizon. 
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The Euclidean action and the entropy behavior of Einstein-dilaton-GB solution was also examined. It was found 
that the inverse temperature is not far from the corresponding Schwarzshild one. When becomes quite large, the 
entropy S is quite equal to the Schwarzshild one Ssw- During decreasing r/i, the difference between S and Saw 
becomes larger. S has always the positive signature. This conclusion does not contradict to Mignemi's formula M, 
obtained by the perturbative method in the case of rather large rh, which has the following form in our variables: 

S = 4irM 2 (l + \-^+\ 2 — -L). (34) 
V M 120 M 4 J V ' 

Our entropy value d32j ) was obtained with the non-perturbative method. So, it is possible to note, that the formula 
(p4|) is indeed correct with the quite large values of and gives the increasing mistake if entropy becomes small. 
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VI. APPENDIX: THE METHOD OF INTEGRATION OVER A PARAMETER OF ORDINARY 
DIFFERENTIAL EQUATIONS GIVEN IN A NON-EVIDENT FORM 

Let us consider a system of usual differential equations with unknown functions Xi(t) S R n which are not solved 
relatively derivatives x^it) to have the following form: 



3=1 
n 

^x' J (t)a 2 j(x,t) = b 2 (x,t), 

3=1 

(35) 

n 

^2xj(t)a n j(x,t) = b n (x,t), 

3=1 

where the coefficients 6j, a,j £ C 1 [R n x R 1 }. For the simplicity it is convenient to introduce (n x n) matrix 

A(x,t) = {cLij(x,t)}ij=i t ... tn . The system ( pq ) is a linear x'^t) system, therefore, it can be solved and its solution has 

the following matrix form: 

x'(t) = yl" 1 (x,t)& t (a;,i), (36) 

where transposed vector b^(x,t) is equal to (bi(x,t), b%(x, t),..., b n (x, t)y. Let to to be the unit value of t and 
x(to) = Xq. So, one has the Cauchy problem that is to continue the solution of ( ^B|) to the manifold t € [to,ti]. 
According to the existence theorem, the solution x(t, £q,xq) of ( ^5[ ) would exist only in the case of the main system 
discriminant to be not equal to zero on the solutions of (|3q), beginning from the point (xo, to), i.e.: 

detA(x(t,t ,x ),t) ^ 0. (37) 



If in any point t* the following condition is realized, i.e. det A(x(t, to, xq), t) = 0, the system ( |3q ) has a non-single 
solution, and a numerical integration performed by the usual Runge-Kutta, Adams, ets. method will stop before the 
point t* . So, it is necessary to choose such a numerical integration strategy which allows to pass the point t*. The 
same method used for other problem classes was discussed early, for example, by V.A.Egorov et all Q. 

For the purposes discussed in the paper one can introduce a new independent variable s and further to consider x 
and t to be functions of s: x — x(s), t — t(s). Therefore, 

§ = ?§=«r 1 - (38) 

at as at 
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Putting (|38|) to the system (|35|), one obtains a new homogeneous system which consists of n equations with (n + 1) 
unknowns. One needs to add (n + l)th equation representing the normalization condition to the system, which takes 
the following form: 



(39) 



3=1 S S 

Edxj . . dt , . 
-^a 2j (x,t) - -^h(x,t) = 0, 

3=1 S S 

Edxj , , dt , . 
—a nj (x,t) - —b n (x,t) = 0, 

3=1 

±^w -§*<•)-»■ 

3=1 

where new functions r(s) € R n+l . These new functions r(s) are fixed by the following normalization conditions: 

, . dxi/ds , . dt/ds , , 

n (s)= 11 n(s)= 1 (40) 

n In 

J E (dxt/ds) 2 + (dt/ds) 2 / £ (dxi/ds) 2 + (dt/ds) 2 

So, r(s) is the tangent vector to the solution trajectory in (n + 1) phase-space, therefore, the integration is proceeded 
along the solution trajectory, not along the radial coordinate r, as usual performed. At a first step the most convenient 
choice of f is the following: Tj — 0, r t — ±1, where its sign depends upon the solution direction moving. At the 
next step r[s) are extrapolated by the Legendre polynoms using the r(s) values from the previous steps according to 
formulae ([Io|). Therefore the system (|39| ) remains the linear one. One can argue that if the integrating step vanishes 
the solution of the system (|39|) w ith the extrapolating functions T(s)reduce to the solution of the system (^) with 
the functions f(s) defined in (ffO|). In such case our system (^J), (|J), (|9|) is written as follows: 
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kg = 0, 



= 0, 



T — = 1 

^ ds 1 ' 



(41) 



where stroke denotes d/dr and , 6j have the form: 
an = 0, 

o u = -m^ 2 r + 4e- 20 A<?! ) '( ( T 2 - 3A), 

o 13 ^4e- 20 ACT(A-o- 2 ), 

a 21 = m^crV + 4e- 2l ^A0'2Acr, 

a 22 = -m 2 ,cr 2 (AV + 2 A) - 4e~ 20 A0'2A3A', 

a 23 = 4e" 20 A2AA'cr, 

a 31 = 4e" 20 A(Afj - a 3 ), 

a 32 = 2m 2 pl cr 2 Ar 2 <j>' + 4e- 2<^i A(-3AA , + AV 2 ), 
a 33 = -2m 2 ,cr 2 Ar 2 <7, 
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h = -m 2 pl a 3 r 2 (<t>'f + 4e- 2 ^Aa(A - a 2 )2(0') 2 , 

b 2 = -m 2 pl a 2 (2A'a + 2Ara{4>') 2 ) + 4e^A2<jAA'2(< ? /) 2 - 4e- 2 ^A0'2(A') V, 
63 = 2to 2 / ct 2 (AV 2 ct< ? !)' + 2Ara(f>') - Ae^ X(A') 2 a. 

The integration is proceeded by the independent variable quantity s and starts from the point of s = 0. At the first 
step it is necessary the following condition to be fulfilled: det B=det A ^ 0, where B is the equation matrix defined 
as follows: 

/ 1 
10 

d _ an ai2 
a 2 i a 2 2 
a 3 i a 32 

At every integration step one calculates f by the formulae (Q) . The 8th order Runge-Kutta program for integration 
by one step was used in our work. This program divides the given step automatic ally to reach the accuracy required. 

One can argue that if the integration step becomes quite small, the matrix B ( f42|) will not be degenerated at that 
trajectory point and rangA > (n — 1). When rangA < (n — 1), a ramifying point for two or more branches can exist 
in that point. Such a point can not be passed by the describing strategy. 

Matrix B can be turned out by the Singular Value Decomposition (SVD) method It decomposes the given 
matrix B to the product of three matrixes UTiV\ where U and V are the unitary matrixes (n+ 1) x (n+ 1) dimension 
and £ is the diagonal matrix of absolute own values of B, i.e. S = {da}, <Ju > 0. Therefore, B^ 1 = V'E~ 1 U' ! where 
= {l/<7jj}. The most effective SVD algorithm can be found in the Numerical Analytical Group (NAG). It is very 
useful because during the calculation one sees the degeneracy's order of matrix B (singular values an vanish). 

In the calculations it is desirable to work with Plank unit values, setting m p i — 1. For simplicity the parameter A 
in the action (|l|) was also set to be equal to one. Equations fl2C|), ( |3l| ) were added to the main system (El]) to reach 
the maximum available accuracy for calculating ADM mass M, dilaton charge D and Euclidean action 
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-h 

-h 
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FIG. 1. The dependence of the metric functions A (a), a (b) and dilaton function exp( — 2<f>) (c) versus the radial coordinate 
r when the event horizon value ru is equal to 10.0. 



FIG. 2. 3D plots of the dependence of the metric function A(r) (a) and the dilaton function exp(— 2<f>(r)) (b) versus the 
event horizon value ru- 



FIG. 3. The entropy behavior of the Schwarzshild solution (line) and GB one (squares) versus the event horizon value rn- 
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TABLE I. The dependence of ADM mass M, dilaton charge D, GB inverse temperature j3 and Schwarzshild one Psw versus 
the event horizon value rn- 



Th 


3.5 


5.0 


10.0 


15.0 


20.0 


30.0 


100.0 


M 


1.8531 


2.5365 


5.0048 


7.5014 


10.0006 


15.0001 


50.0000 


D 


-0.4658 


-0.3632 


-0.1951 


-0.1319 


-0.0994 


-0.0665 


-0.01999 


(3 


43.1675 


62.8607 


125.6630 


188.4898 


251.2456 


376.9911 


1256.6370 


Psw 


43.9824 


62.8320 


125.6640 


188.4960 


251.3280 


376.9920 


1256.6400 
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